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, Abstract. We examined the properties of a sample of BATSE Gamma-Ray Bursts (GRBs) 

, comprising events which have indications of association with a supernova (SN), some on 

■ the basis of indications of re-brightening in the optical afterglow light curve, but in most 

o ■ 

■ cases based only on the 'loose' temporal and directional coincidence inferred from the cross 

correlation of catalogs. Despite of the large uncertainties in the latter selection method, the 

Q ■ temporal and spectral analysis reveal three interesting statistical results when the sample 

i-t ' 

■J-J ' is compared with that of all the BATSE GRBs: the GRBs tentatively associated with SNe 

are found to predominantly (in ~ 80% of the cases) have single-peaked light curves, a 

^ ' softer spectrum (i.e. low energy power law index a 1.5) and tend not to follow the 

' Lag-Luminosity and Isotropic Energy-Peak Energy correlations. These three independent 

' statistical properties point toward the existence of a significant number of under-luminous, 

GRB 980425-like events constituting - at least from an observational point of view - a tail 
or a separate class with respect to the whole of the BATSE GRB events. The unusually high 
percentage of SN Ibc among those identified by the catalog cross-correlation (factor ~ 4 
higher than expected from SN catalog statistics) reinforces the non-randomness of (some 
of) the selected events. 
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1. Introduction 

An increasing number of recent observations support the association of (at least) some GRBs 
with SN events, stimulating both the theoretical understanding/modeling of GRB progenitors 
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and intensive observational campaigns aimed at finding further and more detailed evidence of 
associations. 

GRB 980425 was the first burst with tantalizing evidence for a physical association with 
a SN (SN 1998bw, Galama et al. 1998). The association of the two events has been discussed 
in the frame of models involving the core collapse of massive stars (e.g. Woosley et al. 1999; 
MacFadyen & Woosley 1999; Vietri & Stella 1998). However, the unusual properties of both 
events within their classes - i.e. the high luminosity of the SN, the unusually low luminosity 
and variability of the GRB y-ray emission and the lack of detected X-ray afterglow - as well as 
the uniqueness of this association, have not allowed to draw firm conclusions about the possible 
progenitor for a few years. 

More recently two further GRB/SN associations based on the spectroscopic detection of a SN 
spectrum in the GRB afterglow added robust evidence to the association of the two phenomena. 
GRB 030329 (Stanek et al. 2003; Hjorth et al. 2003) has been associated to SN 2003dh from the 
appearance of typical SN Ic spectral features in its optical spectrum and the (under-energetic) 
GRB 03 1203 revealed an optical spectrum consistent with the contribution of an underlying type 
Ic SN (SN 20031w; Malesani et al. 2004). 

Despite of these promising findings, the paucity of cases of GRB directly associated with SN 
Ib/c, some peculiar properties of both the GRB and SN events in these few robust associations 
and the uncertainties on the event rates in both classes (see however Podsiadlowski et al. 2004) 
do not allow a clear insight and different possibilities are still open: selection effects might play 
a role in detecting SN features, there might be an (inverse) relation between the energetics of the 
two events (under-luminous GRB vs Hypernova event), or bursts associated with core collapse 
SN could represent a subset of the whole population, which might be even characterized by 
properties (e.g. energetics, duration, temporal variability) intermediate between the two class of 
events (e.g. Delia Valle 2004 for an overview; Soderberg 2004). 

Given the limited number of 'direct' associations (i.e. with evidence of SN spectral features 
in the GRB afterglow), in order to test for the nature of any physical connection and thus in- 
vestigate the progenitor(s) of GRB, other more indirect indications could be considered to select 
candidates GRB/SN associations. These include the re-brightening in the optical afterglow light 
curve decay, possible signature of SN emission emerging at times when the GRB has sufficiently 
faded (Bloom et al. 1999; Germany et al. 2000). Indeed about ten such associations have been 
claimed on such a basis (Zeh et al. 2004). Interestingly Zeh et al. (2004) found indication of 
re-brightening in all GRBs with known redshift z < 0.7. Similar evidence (i.e. a bump in the 
late time light curve) have been found also in the case of a few X-ray Flashes (see Fynbo et al. 
2004, Soderberg et al. 2005). An even more indirect and uncertain, but possibly more efficient 
approach, consists of systematic searches for directional and temporal coincidences of GRBs and 
SN from the cross correlation of GRBs and SN catalogs (Wang & Wheeler 1998, Hudec et al. 
1999). These searches are clearly affected by i) the large uncertainties in the GRB position (typ- 
ically of a few degrees) (e.g. Kippen et al. 1998) leading to possible serious overestimates of the 
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number of real associations, but also ii) the completeness and flux limit of the SN/GRB catalogs 
used. Nevertheless, if the indications by Zeh et al. (2004) prove to hold, catalog cross-correlations 
should indeed find evidence for a significant number of real associations. 

In principle a more powerful method for selecting candidate associations from the catalog 
cross correlation could be envisaged if GRB associated to SN were characterized by particular 
temporal and/or spectral properties with respect to the whole of the GRB population. Indeed 
Norris et al. (1999) explored the hypothesis that there exists a subclass of GRBs related to SN 
whose y-ray prompt light curve presents a single smooth peak, similar to that of GRB 980425 
(as suggested by Bloom et al. 1998), but with negative results. 

Here we focus on examining the possible peculiarities of the candidate events by investi- 
gating the temporal and spectral properties of GRBs for which there is an indication of being 
associated with a SN as a "class". More specifically we considered direct and indirect (tentative) 
associations, the latter ones based either on the light curve re-brightening or on the positional- 
temporal coincidence from catalogs cross correlations. To this aim we searched the literature 
for the most up-to-date sample of GRB/SN candidates including the three robust cases of asso- 
ciation (GRB 980425, GRB 030329, GRB 031203) and the 11 cases where the afterglow light 
curve re-brightening gives evidence for a coincident supernova explosion, that is GRB 021211 
(Delia Valle et al. 2003), GRB 011121 (Gamavich et al. 2003), GRB 020410 (Levan et al. 2004), 
GRB 041006 (Stanek et al. 2005) and 7 GRBs examined by Zeh et al. (2004) (GRB 970228, 
GRB 980703, GRB 990712, GRB 991208, GRB 000911, GRB 010921, GRB 020405). From 
the available data we find that the light curve of these events consists of a single peak in 50% 
of cases and that their spectral parameters indicate a soft low energy spectrum. Given these sug- 
gestions we then extended the number of candidates by including those from both past catalogue 
cross-correlations and updating them using more recent GRB and SN catalogs. 

Clearly, if GRBs in the sample turned out to be statistically (although not individually) char- 
acterized by any peculiar property - with respect to a random GRB sample - this would have 
interesting implications: i) it would indicate the existence of a SN-related sub-class of GRBs; ii) 
it would provide an empirical tool for efficiently identifying such events; iii) most importantly, it 
would offer clues on the physical understanding of their origin (indicating a mechanism different 
from that underlying GRB without SN signatures). 

The outUne of the paper is the following. The selected GRB sample is described in Section 2 
and its general properties are discussed in Section 3. In Sections 4 and 5 we present the detailed 
temporal and spectral analysis. We discuss in Section 6 the reliability of the GRB/SN associations 
and the nature of the selected GRB by means of (two of) the phenomenological intrinsic rela- 
tions found for GRB with estimated redshift. Section 7 summarizes our conclusions. We adopt a 
ACDM cosmology with Qa =0.7, Qm =0.3, //q = 65 km s"^ 
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2. Sample selection 



Let us describe in detail the (heterogeneous) sample considered here, which basically includes 
the BATSE GRBs with indications of association with SN and for which temporal and spec- 
tral information could be retrieved, for a total of 36 GRB/SN. The whole sample with its basic 
properties is hsted in Table 1 and for simphcity hereafter referred to as the GRB/SN sample. The 
three robust associations based on the spectroscopic identification of a SN, namely GRB 980425, 
GRB 031203 and GRB 030329 will be considered as reference cases throughout this work. 

More precisely the GRB/SN sample comprises: 

i) 7 GRB/SN suggested by a bump in the late-time optical afterglow (Bloom et al. 1999, Zeh et 
al. 2004), robust spatial and temporal association (GRB 980425) or an indication of a line feature 
in X-rays (Piro et al. 1999), or exceptional SN properties (high luminosity, peculiar spectrum; 
e.g. Germany et al. 2000). Of all the claimed cases we included those with available BATSE data; 

ii) 15 GRB/SN from the previous searches of temporal and spatial coincidence based on the 
cross-correlation of BATSE and SN catalogs: 11 pairs by Wang & Wheeler (1998) and 4 pairs 
byHudecetal. (1999); 

iii) 14 GRB/SN identified in this work from an update of the catalog cross-correlation (including 
approximately the dates previously not covered, after Nov. 20, 1997 up to May 26, 2000). In this 
period the BATSE GRB catalog ^ contains 721 GRB and the SN catalog of lAUC^ reports 495 
SN. 

The catalog cross correlation method that we adopted is based on the following criteria. 
For the positional coincidence we required a SN to be within the BATSE GRB error box to 
which the constant 1.6 degrees of r.m.s. systematic error has been added (see e.g. Meegan et al. 
1996). The exact SN phase at detection is not known (typically the optical maximum is reached 
in ~20 days) and furthermore no robust guess is possible for the delay between the SN and 
GRB explosions, as any estimate is strongly model dependent: ~ a few months up to few years 
are allowed in the supranova model (Vietri & Stella 1998) while the explosions occur almost 
simultaneously in the coUapsar model (e.g. MacFadyen & Woosley 1998). Although the latter 
case is theoretically and observationally favored (e.g. Delia Valle 2004 for a review), a large 
temporal window has been allowed for the possible correlation of the two events, requiring the 
SN date to be between -10 to -1-30 days with respect to the GRB trigger. No further constraints 
were imposed in the selection of either events and thus the results are biased by the SN sky 
distribution and the detection incompleteness at faint magnitudes; we considered both long and 
short GRB and no limit was imposed on the GRB count rate. Also, despite of the indications 
from the three robust associations, we did not impose any selection on the SN type in our search 
(see however Section 5). 

' http://cossc.gsfc.nasa.gov/batse/BATSE_Ctlg/index.html 
^ http://cfa-www.harvard.edu/iau/lists/SN.html 
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Note that in the previous catalogs cross correlations analysis different search conditions were 
imposed. Wang & Wheeler (1998) searched for SN with coordinates within 2cr BATSE error 
boxes and supposed different time delays (according to the SN type) between SN and GRB. They 
also limited their search to SN at z < 0. 1 . As this limit, however, would exclude 2 (GRB 030329 
and GRB 031203) out of the 3 more robust associations we did not impose any redshift cut. 
Hudec et al. (1999) reported the correlations for SN detected within the GRB error box up to 30 
days after the GRB. Thus the compiled sample is not homogeneous as different criteria were ap- 
plied. We considered re-examining the associations for the periods examined by Wang & Wheeler 
(1998) and Hudec et al. (1999). However the BATSE instrumental trigger parameters have been 
changed during the BATSE 9 yr activity (Paciesas et al. 1999), until February 1997, when they 
have been fixed. Therefore, a re-analysis would have been affected by this further bias, and thus 
not led to any increased completeness of the GRB/SN sample. 

This updated cross-correlation yields 16 new GRB/SN coincident events. However, coordi- 
nates from the IPN^ for two GRB revealed that the associations were spurious (and thus excluded 
from the sample). In Fig. 1 error boxes and time delays are displayed for 31 GRB/SN (2 GRBs 
had multiple associations and 3 GRBs had an indication of an underlying SN from the afterglow 
observations). 

The significance of such blind (catalog based) searches was studied by Kippen et al. (1998), 
who assigned a significance to the observed correlation statistics by generating a set of random 
GRB catalogs. They find no evidence for association for GRBs having precise BATSE/Ulysses 
localizations, although they could not exclude it for a subset of weak bursts. While these results 
point toward a low chance to select real coincidences with such a procedure, the peculiarity of 
some properties of the GRB/SN sample (shown below) is instead suggesting that such method 
does not select a random sub-sample of BATSE GRBs. 

3. GRB/SN sample properties 

Based on the properties of GRB 980425 it has been suggested (see e.g. Bloom et al. 1998) that, if 
there is a physical connection between a sub-sample of GRBs and SN, such GRBs might exhibit 
distinct characteristics. Therefore we have compared the global properties of the GRB/SN sample 
with those of the entire BATSE sample. 

No qualitative (nor statistical) differences are found as far as duration (represented by Tgo), 
spectral hardness (HR32) ratio and total burst fluence are concerned (see Fig. 2, top and bottom 
panels), as both the range and shape of the GRB/SN distributions are indistinguishable from 
those of the rest of GRBs in the BATSE catalog. 

As expected from the typical worsening of the positional accuracy with decreasing photon 
statistics (see Kippen et al. 1998), the GRBs identified by the catalogs cross-correlations are 
mostly weak, with typical peak (50-300 keV) flux < 1.5 phot cm"^ s"^ (see Fig. 2). 
^ http://www.ssl.berkeley.edu/ipn3/interpla.litml 
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Fig. 1. BATSE error boxes and time delays for the GRB/SN sample. The time delay is between 
the SN discovery date and the GRB event. 

While no peculiarity is found for the above quantities in the GRB/SN sample, in the next 
Section we present a more detailed analysis of their GRB temporal profile and spectra'*. Notice 
that the bias on the flux as well as the selection criteria of the GRB/SN sample result in a "non 
optimal" statistics for the temporal and spectral analysis. This in turn requires a careful selection 
of a BATSE GRB "control sample", as clearly previous analysis of BATSE GRBs have been 
typically performed on brighter samples. 

4. Light curves analysis 

In this Section we present the analysis of the temporal profiles of the GRB/SN sample and com- 
pare the results with the (known) properties of larger GRB control samples. 

While the temporal profiles of GRBs display a wide diversity (see e.g. Norris et al. 1999) 
we recall for future reference some of the important systematic findings on their properties: (1) 
flux peaks occur later at lower energies (spectral lags - Norris et al. 1996, Norris et al. 2000); 
(2) the pulse parameter distributions are consistent with lognormal distributions (Quilligan et al. 
2002); (3) less variable GRBs are found to be intrinsically less luminous (variability-luminosity 
correlation, Fenimore & Ramirez-Ruiz 2000). 

4.1. Analysis 

For most of the GRBs (16/36) we analyzed the time profiles from the BATSE concatenated 64- 
ms data derived from the eight BATSE Large Area Detectors (LADs); these data represent the 
detector count rates in 4 energy channels (25-50, 50-100, 100-300 and >300 keV). We fitted both 
the energy-integrated flux and the individual 4 energy channel light curves, according to the 

* Note that in the following the GRB/SN events reported in each case/figure is different, depending on 
the quality of the available data. 



1 A 

A 

A a 



Z. Bosnjak et al.: GRBs associated with SNe 



7 



1000.00 r 



100.00 r 



m 10.00 



1.00 r 



0.10 r 



1 •♦• 
.. • » 



0.01 I , , , 

0.1 1.0 10.0 100.0 

peak flux / phot cm"^ s"' 



> 

O 

o 



100.0 r 



10.0 r 



> 

O 

c 



1.0 r 



c 
o 



0.1 L 

10" 



7 6 5 4 

St.i(25-* keV) / ergs cm"" 



10" 



Fig. 2. Comparison of the durations Tgo, peak fluxes, hardness ratios (i.e. HR32) and total flu- 
ences of the GRB/SN and the entire BATSE samples. Triangles, squares and circles indicate 
the tentative associations found by Wang & Wheeler (1998), Hudec et al. (1999) and this work, 
respectively. Diamonds are the GRB with an indication of a SN event in the afterglow. FiUed 
symbols represent long GRB {T^o >2 s). 



count statistics available for each event. The duration of 6 events was shorter than 2 s: in these 
cases we could thus analyze the Time-Tagged Event (TTE) data with 2-ms time resolution (which 
are bufl'er-size Umited^) integrated over the 4 energy channels. In the analysis of these data type 
we also accounted for the «30 ms ofi'set between trigger time (as reported in the BATSE catalog) 
and TTE data load time. For the remaining 14 GRB/SN events either no reprocessed light curves 
were found in the pubUc archive (9 cases) or the low statistics of the Ught curve did not allow 
any reasonable fit (5 cases). 

First we determined for every GRB the background count rate selecting two time intervals 
not containing the burst^, and modeled it with a 2nd order polynomial. We then built a pulse 

^ http://cossc.gsfc.nasa.gov/batse/batseburst/tte/index.html 

* http://cossc.gsfc.nasa.gov/batse/batseburst/sixtyfour_ms/ bckgnd_fits.html 
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model by adding the best fit background model to the lognormal pulse profile (Brock et al. 1994) 
and fit the light curve as: 







1 / iog(t-to)-^ y 



, t>tQ 

Otherwise 



where fi and cr are the mean and standard deviation, respectively, and the time is scaled to the 
pulse onset to ^. Two reasons are behind this choice of profile: 1) most of the GRBs in the sample 
appear as single-peaked (on the 64 ms and 2 ms timescale for long and short bursts, respectively); 
2) it allows a direct comparison with the results of previous studies (e.g. Norris et al. 1999). The 
best fits (x^ ~ 0.83 - 1.18) are obtained for the integrated light curves, and in 11 (out of total 22) 
cases an acceptable fit was obtained for individual channel light curves. 



4.2. Results 

The resulting distribution of cr is shown in Fig. 3 (top panel) for all the pulses in the light curves 
of the 22 GRBs for which the fit was performed. The distribution of cr vs. /i is reported in the 
bottom panel. As expected short bursts are characterized by cr lower than long ones and their 
pulses are less asymmetric. While the statistics is rather poor, we note a global tendency for an 
anti-correlation between cr and yU both in long and short GRB/SN: there is no pulse with large 
(i.e. long rise time) and large cr (i.e. very asymmetric), namely they are typically FRED-like (i.e. 
Fast Rise Exponential Decay). 

In order to compare the GRB/SN pulse properties with a control sample we examined also 
the distribution of their FWHM (64 ms time resolution). In Fig. 4 this is shown together with 
the distributions derived by Norris et al. (1999) and Quilligan et al. (2002). The latter, referring 
to the brightest (long) 319 BATSE GRBs, is very broad, peaking around 0.7 s, and significantly 
different from that of the GRB/SN sample. Interestingly, however, the FWHM distribution for 
the subset of long GRB/SN events (solid line in Fig. 4) appears to be similar to that derived by 
Norris et al. (1999) from the analysis of a sample of long, single -peaked GRB. Indeed a K-S 
test returns (for long, single -peaked GRB/SN) a probabiUty Pj5:s=0.98 that the two data sets 
originate from the same parent distribution. 

In Fig. 5 (top panel) the pulse amplitudes (only for 64 ms resolution timescale) vs. FWHM 
are displayed. Quilligan et al. (2002) find a weak negative correlation between these parame- 
ters. While there is no evidence for a trend for the whole GRB/SN sample, an anti-correlation 
(Spearman rank correlation coefficient of -0.65 and probability 0.003) is present among the 
GRB/SN selected on the basis of the cross-correlation only. The interpretation of this is not 
clear: it might be simply due to the bias toward low fluxes in GRB/SN selected from the cata- 
^ The lognormal pulse model is skewed to the right (as determined by cr). 

* While formally the results by Norris et al (1999) refer to fits of count rates from channels CH2-I-3, 
typically the signals in CHI and CH4 are weaker, thus not significantly affecting the comparison. 
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Fig. 3. Top: Distribution of cr for the lognormal pulse model. The solid and dashed lines represent 
the parameters for long bursts and the contribution of short bursts. All the pulses in the time 
history are represented. Bottom: cr vs fi, i.e. the mean in log time. Symbols as in Fig. 2. 
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Fig. 4. Distributions of full widths at half maximum. The solid/dashed lines represent the contri- 
bution of long/short GRBs. The dotted line refers to the single-peaked bursts analyzed by Norris 
et al (1999), while the dash-dotted line are the bright, long bursts analyzed by Quilligan et al 
(2002). For clarity the distributions are normalized to the maximum number of data points (i.e. 
long bursts of Quilligan et al 2002). 
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Fig. 5. Top: Amplitude (fitted at 64 ms resolution) vs. FWHM as obtained from the fit to the 
total counts in the four channels. Bottom: Time lag between fitted peaks (in CH2 and CH3) vs. 
FWHM. Only long bursts are shown. Symbols as in Fig. 2. 

log cross-correlation and/or indicate that indeed the latter are largely spurious associations, not 
representing (atypical) properties of GRB physically associated with SN events. 

A tendency for wider pulses to have larger time lags appears instead for the GRB/SN sam- 
ple (Fig. 5, bottom panel). There the spectral lag (i.e. the delay between the peaks of the low 
with respect to the high energy channel) is plotted against the FWHM of the fitted pulse. This 
behavior is consistent with what reported by Norris et al. (1996) on the basis of a sample of 46 
bright, long bursts: they found that the lag results mostly from the shift in the pulse centroid and 
thus the correlation is induced by the narrowing of the pulses toward higher energies. Different 
interpretations have been proposed for such property (e.g. Salmonson 2000, Schaefer 2004) and 
very recently it has been suggested that it might be accounted for in the jet-viewing-angle sce- 
nario (loka & Nakamura 2001) where both (small) pulse amplitudes and (long) spectral lags are 
related to (large) viewing angles of the jet. As shown in the top and bottom panels of Fig. 5, in the 
case of the GRB/SN sample (but excluding the 7 GRB/SN determined by the light curve bumps) 
small ampUtude pulses could correspond to large spectral lags. 

To summarize: within the limited statistics and differences in the selection criteria of the 
GRB/SN sample, the GRB/SN events appear to be broadly similar to the whole of the BATSE 
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GRB population. However, the analysis of the FWHM distributions highlighted one intriguing 
possibility, namely that the GRB/SN events behave more typically like single-peaked GRBs. 
Although aware that the bias toward weak GRB in the sample might be responsible for the above 
finding - Norris et al. (1999) already cautioned that dim bursts might have a multiple pulse 
structure which is difficult to detect - we explored this aspect in detail. 

4.3. A connection between single-peaked GRBs and SN? 

As already mentioned, Bloom et al. (1998) suggested that GRBs associated with SN might have 
distinctive temporal and spectral properties, similar to those of GRB 980425, namely a smooth, 
single peaked, long-duration (790=34.8 s) light curve and faint emission above 300 keV. Norris et 
al. (1999) considered this issue and searched the BATSE catalog with these criteria, but no other 
event similar to GRB 980425 and spatially/temporally correlated with any known SN was found. 
In our GRB/SN sample we instead found that the large majority (~ 80 %) of GRB are single 
peaked. Of those not selected by the catalog cross-correlation 5 out of 7 GRBs are single-peaked 
events. It is worth stressing that of the other two spectral associations, GRB 031203 shows a 
single peaked light curve (Soderberg et al. 2004), while GRB 030329 presents different temporal 
properties, i.e. a profile with two peaks of similar intensity (Vanderspek et al. 2004). 

In order to quantitatively test whether the GRB/SN sample comprises an unusually high 
number of single-peaked GRBs (1) it has been necessary to determine an objective way of char- 
acterizing the "single-peakness" and (2) due to the limited size of the GRB/SN sample (29 GRBs 
with 64 ms data) the probability of finding a similar large percentage of single-peaked events 
had to be computed via random extractions of GRB sub-samples from the BATSE catalog. 

4.4. Method 

We thus built an algorithm that identifies single-peaked bursts in the BATSE catalog. As the time 
profiles of GRBs difl'er to great extent it is difiicult to have a unique approach applicable to every 
burst: eventually we imposed conditions just to determine whether the burst has one or more than 
one peak. 

Li & Fenimore (1996) proposed a method for identifying peaks in GRB time histories, and 
a modification of this was also applied by Nakar & Piran (2002) in the temporal analysis of 
short bursts. It consists of defining a count bin Cp at time tp as a candidate peak if it satisfied the 
condition: 



where Ci 2 are neighboring count bins at times ti^2- Ny^r determines how rigorous the search is: 
too large A^var might clearly result in an algorithm insensitive to faint peaks; too low A'var could 
identify statistical variations as peaks. The previous works showed that the choice of Nyai=5 
optimizes the algorithm. 




(1) 
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Fig. 6. Distributions of the number of single-peaked GRB in sub-samples of N-29 GRB. The 
dotted, dashed and solid histograms represent the distributions for respectively n-5Q, 150, 250 
extractions out of the total sample of 1445 bursts. The lognormal fit corresponds to the case 
n=25Q. 

The search was performed on the background-subtracted light curves, and thus a linear back- 
ground was fitted to data before and after the burst trigger, and subtracted. The resulting time 
histories were biimed to 128 ms resolution. 

As one of the problems encountered was that sharp structures in the time history could be 
mis-identified as peaks, we applied a wavelet analysis for determining the dominant modes of 
variability. Wavelet transform is convenient to use when the analyzed signal consists of short 
spikes, that are not well approximated when doing a Fourier analysis: the wavelet coefiicients 
resulting from the convolution of the signal with the wavelet basis functions represent the signal 
variation at a particular resolution (see e.g. Kolaczyk 1997, QuiUigan et al. 2002 for details). We 
reconstructed the original signal by keeping only the maxima of the wavelet transforms; a very 
limited number (~ 25) of the wavelet coeflftcients was used to reproduce the signal. Although this 
does not guarantee the generation of a smooth light curve, with this approach the unique peak 
finding condition could be applied for the whole set of GRBs. 

The search for peaks was limited to the burst region, excluding periods before the burst which 
may contain a precursor. We tested several choices for the condition of peak finding, optimizing 
it to recognize what would be visually identified as a single-peaked GRB, and implemented 
the slightly modified peak-finding conditions: a) adopted the value A'var=7, and b) added the 
condition for a peak to have at least 20 % more flux than the neighboring 'vaUeys' . 

4.5. Results 

We applied the algorithm to the GRB/SN sample and this defined 23 (out of 29) single-peak 
bursts: this difl'ers by less than 10% from the number we find by visual inspection. 
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To establish the chance probability of such occurrence, we selected GRBs from the BATSE 
catalog having fluences in the range span by the GRB/SN events (between 6.9 x 10"** and 
8.4 X 10"* erg cm"^), ending up with 1587 GRBs, and systematically applied the algorithm 
to their 64-ms resolution data (for 142 out of them the algorithm did not identify any peak, as 
a consequence of the rigorous conditions imposed). We then randomly extracted sub-samples of 
the same size (N-29) of the GRB/SN sample, repeating the extraction n times. In Fig. 6 we show 
the distribution of the fraction of single-peak GRBs for n-50, 150 and 250. 

We fitted a lognormal function to the probability distribution^, i.e: 

Pix) = _}_e^-(Mx)-Mf/2cr^)^ (2) 

for the case n =250 and find a 0.4 % probability of having 23/29 single-peaked bursts. If the 
number of single-peaked bursts in the sample were 10% lower (i.e. 20 or more events) the prob- 
ability of observing it would stiU be ~ 5%. Note that the distribution has a maximum for ~ 15 
single-peaked bursts; this is consistent with the result by Norris et al. (1999) who found 68/116 
bursts with at least one pulse - suggesting however that it is difficult to assign an exact number of 
pulses as they observed a continuum of lower emissions superposed on otherwise single-peaked 
events. 

The homogeneous treatment of the GRB/SN and random control sub-samples should ensure 
that the above result is unbiased against effects that may vary the number of detected pulses in the 
GRB time history, such as the rigorous conditions on peak candidates, the fluence range, the (still 
present) spikes mis-interpreted as pulses and the presence of short GRBs which in fact appear to 
be multi-peaked when seen at 2 ms resolution (there was 23% of short bursts in the sample). 

5. Spectral analysis 

Let us now consider whether GRB/SN show any peculiarity in terms of spectral properties. We 
recall that the spectrum of GRB 980425 is 'unusual', peaking at ~ 100 keV (Jimenez et al. 2001), 
below the typical 300 keV peak energy. Also GRB 030329 is characterized by a similar low peak 
energy (Vanderspeck et al. 2004), while GRB 031203 has a larger peak energy (> 212 keV, 
Sazonov et al. 2004). 

GRB spectra in the BATSE energy range (25-1800 keV) are usually described via the phe- 
nomenological Band representation (Band et al. 1993), which consists of two smoothly con- 
nected low- and high- energy power laws. Despite of a non-negligible dispersion of the fitted 

parameters, the typical values for the low and high energy power-laws photon indices are a 1 

and /3 2.5, respectively, and the spectral peak energy fipeak (i-e. the peak in a vFv representa- 
tion) is around a few hundred keV (Preece et al. 2000). The Band model in general successfully 
represents both the time integrated as well as the time resolved spectra (Ford et al. 1995). 

' Lognormal distributions have been already used to describe the distributions of various GRB properties 
(see e.g. Li & Fenimore 1996, QuiUigan et al. 2002) 
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Following the widely adopted methodology described by Preece et al. (2000), we analyzed 
the GRB/SN spectra obtained by integrating the signal over the burst duration (a time resolved 
spectral analysis was possible only in a few cases due to the low signal-to-noise ratio). 

5.1. Analysis 

The spectral analysis was heavily limited because in many cases the GRB signal was too low 
for any reasonable spectral fit. For just over half of the bursts (16) we used the High Energy 
Resolution Burst (HERB) data from the most brightly illuminated BATSE LAD (128 energy 
chaimels). For short bursts or in those cases with no good time coverage with the HERB data 
we analyzed the Medium Energy Resolution (MER) data (5 cases), which consist of 16-energy 
chaimel spectra with a time resolution of 16 ms during the first half of the total duration and 64 
ms during the second half. 

The background spectrum was extracted and averaged from two time intervals, before and 
after the burst trigger. The time interval for the spectrum accumulation typically starts at the 
trigger and lasts for the duration as estimated by Tgo. The spectrum was eventually re-binned in 
energy in order to have enough statistics to apply the;^^ goodness-of-fit test 

The spectra were analyzed using the XSPEC v. 11. 1.0 package, performing the fitting with 
the Band model in the ~30 — 1700 keV energy range. In some GRBs this model fit resulted in 
unconstrained spectral parameters and/or unacceptable reduced (18 cases): if mainly due to 
the low statistics at high energies we tested as alternative models the cutoff-power law (CPL; for 
1 GRB) and broken power-law (BPL; 6 GRBs) spectral functions, while if due to a low signal 
in most chaimels, we adopted the simpler single power-law model fit (11 GRBs). The same set 
of models was applied to larger samples of bright bursts (Preece et al. 2000): we stress that 
our aim is to characterize the spectra of the GRB/SN sample rather than establish which model 
preferentially fits them (in fact the choice of a particular model may not be unique even for bright 
bursts, e.g. Ghirlanda et al. 2002). We considered the model fit acceptable when a reasonable 
(ranging from 0.8 to 1.3) and non systematic residuals were obtained. AU the models (except 
for the single power-law) are characterized by a low energy power-law oc E", a high energy 
power-law (Band and BPL models) oc and a characteristic energy Epeak- 

5.2. Results 

The distributions of the spectral parameters are shown in Fig. 7 for the GRB/SN and the Preece 
etal. (2000)'s sample 

We binned in order to obtain a minimum of 5 counts/bin. This (minimal) choice has been tested by 

re-binning the data for a minimum of 10 and 20 counts/bin (e.g. Bevington & Robinson 1992) and the 

corresponding fits did not yield either significantly different nor significantly different fit parameters. 
For the comparison the results of the time-resolved spectral analysis by Preece et al. were averaged, as 

for most GRB/SN it was possible to analyze only time integrated spectra. 
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Fig. 7. Spectral parameter distributions obtained for the GRB/SN (dashed Unes) and GRBs ex- 
amined by Preece et al. (2000) (dotted lines). For clarity the number of GRB/SN events is nor- 
malized to the total number of GRBs in the larger control sample by Preece et al. (2000). In the 
top panel the results obtained by fitting a single power law model are indicated separately (solid 
line). A high energy power-law index yS = -4 refers to the exponential function fit. Errors for 
the peak energy range from 20 to ~ a few hundred keV. The arrows in each panel indicate the 
values (or limits) corresponding to GRB 980425, GRB 030329 and GRB 031203. Note that the 
highest values in the three distributions refer to GRB 990810 and GRB 971221, found by catalog 
cross-correlation. 



For the GRB/SN sample the average spectral index is a ~-1.5, with almost no bursts with 
-1 < a; < 0, which are instead found in the comparison sample. This suggests that on average 
GRB/SN events have low energy spectra softer than typical GRBs. Note that GRB 980425 and 
GRB 030329 have spectra fitted with a —1.2 (Vanderspek et al. 2004) and GRB 031203 is 
modeled with a single power law with a — 1.6 (Sazonov et al. 2004). It is important to note that 
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the single power law model (fitted to ~ half of GRBs) in general gives a value of a lower than 
that inferred from 'convex' models (Band, CPL or BPL). A K-S test performed on the a values 
obtained by fitting the 'convex' models only yields that the GRB/SN sample is significantly 
different from that reported by Preece et al. (2000), with Pks =0.03. 

The distributions of Speak (for the 10 GRB/SN events that had spectrum fitted by Band, CPL 
or BPL) are reported in the middle panel of Fig. 7. The typical peak energy for GRB/SN is ~ 
220 keV, lower than that of bright BATSE GRBs (from Preece et al. 2000) with Pks ^0.04. As 
mentioned, 2 out of 3 GRBs with spectroscopic evidence of association he in the lowest bin of 
this distribution (see Fig. 7). 

Finally we report in Fig. 7 (bottom panel) the distributions of j8. Within the large errors (fi is 
usually not well constrained because of the S/N decrease in the high energy channels), there is 
no evidence of deviations with respect to bright long bursts. 

6. Redshift estimates 

Clearly a crucial test on the physical connection of our GRB/SN pairs would be the determination 
of the GRB redshifts. Up to date, redshifts are known only for ~ 42 GRBs via the optical emission 
lines of their host galaxies or absorption lines in the optical transient spectra. Due to the high 
y-ray flux, unaffected by extinction, and to the limited number of afterglow detections, several 
attempts have been made in order to find a reliable method for deriving distances from the prompt 
y-ray emission only. Unfortunately the large diversity in the temporal and spectral properties of 
GRBs and the related uncertainties in turn result in significant uncertainties in any such method. 
For our purposes however even an indication of the GRBs distance could rule out or reinforce an 
association. 

Several correlations have been proposed to infer the GRB energetics on the basis of observ- 
able temporal and spectral y-ray properties. Fenimore & Ramirez-Ruiz (2000) (see also Reichart 
et al. 2001) found a correlation between the luminosity and a variability indicator for a few GRBs 
with measured z- smoother bursts appear to be intrinsically less luminous. Another distance in- 
dicator was identified by Norris et al. (2000) as a correlation between GRB luminosities and 
spectral lags. These correlations have then been used to estimate redshifts for a large sample 
of GRBs (Fenimore & Ramirez-Ruiz 2000). Another interesting correlation has been found by 
Amati et al. (2002) (see also Lloyd-Rorming & Ramirez-Ruiz 2002) between the intrinsic peak 
energy Epeak and the isotropic total radiated energy Erad- Alternative methods to determine z have 
been proposed: Atteia (2003) identified a method by combining the spectral parameters and GRB 
duration; Bagoly et al. (2003) suggested a correlation between the peak flux ratios and z. 

It is worth noticing that the calibration of all these relations is based on GRBs with known z, 
but always excluding GRB 980425 due to its peculiar energetics and temporal/spectral properties. 
Also GRB 031203 turned out to be an outsider of some of those correlations. Caution should be 
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Fig. 8. Top panel: Relation by Norris et al. (2000). Bottom panel: Relation by Amati et al. (2000). 
In both panels, the GRB/SN events are marked (diamonds) at the values corresponding to the 
redshift of the associated SN, while the trajectories represent the values corresponding to different 
redshifts, in the range 0.01-4. The dashed lines indicate the correlations as determined by Norris 
et al. (2000) and Amati et al. (2002). GRB/SN events which had only upper or lower limits on 
their peak energy are shown only for the values corresponding to the redshift of the coincident 
SN as arrows. 



thus considered in using the above results for candidate GRB/SN events, which might indeed be 
most similar to GRB 980425 and GRB 031203. 

In the following we focus on two of such methods, for which we could infer more reliably 
the corresponding GRB parameters, namely the Lag-Luminosity correlation (Norris et al. 2000) 
and the Peak Energy-Energy correlation (Amati et al. 2002). 

6. 1. Lag-Luminosity correlation 

Norris et al. (2000) found that the spectral time lag t is anti-correlated with the isotropic peak 
luminosity Liso, namely Liso,53 1.3 x (r/O.Ol s)"'^^ for a sample of 6 bursts with known z. In 
order to derive the lag measure for the GRB/SN the same methodology was applied. In particular: 
i) we calculated the cross correlation functions (CCFs) of the time histories in two different 
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energy channels (see also Band 1997), i.e. CHI and CH3, adding tlie apodization intervals before 
and after the burst (apodization intervals had a constant count number, computed as average value 
in adjacent time bins); ii) the analysis was limited to light curve regions with count rates larger 
than 10% of the peak luminosity; iii) the time series were binned in order to have higher S/N (64 
ms data were binned to 128 ms; for short GRBs TTE data were binned to 4 ms); iv) we fitted a 
linear background to every channel separately in the non-burst portions of the data and subtracted 
it from the total counts; v) finally the CCFs were fitted with a third order polynomial around their 
peaks in order to derive an estimate of the time lag, then corrected by a factor l/(l+z). We 
derived the peak luminosity in the 50-300 keV energy range directly by integration of the GRB 
fitted spectrum for each possible redshift. 

The results are presented in Fig. 8 (top panel). Unfortunately it was possible to estimate time 
lags for only 3 bursts (16 GRB/SN have both spectral data and four-channels time histories avail- 
able and of those most have negligible signal in CHI; also in 3 cases the obtained time lags were 
negative and thus not further considered). The curves represent the intrinsic (i.e. corrected for z) 
Liso and t assuming a range z = 0.01 - 4 for the possible redshift. The points on them indicate 
the values of Liso-t corresponding to the redshift of the putative associated SN, when known. It 
appears that the GRB/SN events could be consistent with the Lag-Luminosity correlation only if 
they were placed at high z (>1.5) indicating that the association are spurious. Alternatively, these 
GRBs might not obey such correlation and thus could still be located at the SN redshift. The 
latter (intriguing) possibility does in fact get some credit from the fact that these GRBs would 
then mimics the characteristics of 2 out of the 3 spectroscopic GRB/SN (i.e. GRB 980425 and 
GRB 031203), possibly forming a separate class. 

6.2. Peak Energy-Energy correlation 

The Peak Energy-Isotropic Energy correlation (Amati et al. 2002) has the form Epeak ^raT- 
Similarly to the procedure described above we computed for the GRB/SN the total bolometric 
fluence (see Amati et al. 2002) by integrating (I-IO'' keV) the best fit spectrum (11 GRB/SN 
events). The results are presented in the Speak-^rad plane in Fig. 8 (bottom panel). The position 
of the symbols again indicates the corresponding values for the SN redshift, while the curved 
lines show the values for the interval z=0.01-4. The GRB/SN events for which only lower/upper 
limits of the peak energy could be estimated are shown only at the z of the possibly associated 
SN. We also mark the values for the 3 spectroscopic associations. Only two of the GRB/SN 
could be consistent with the Amati et al.'s relation: GRB 980703 (at z=0.967, Djorgovski et al. 
1998), and GRB 971221 which might be associated with a SN at z « 0.58 (see Table 1). The 
properties of the rest of GRB/SN are not consistent with those predicted by such relation, but 
most importantly in 10 cases (5 of which with well determined parameters) could not obey the 

Note that both time dilation and channel energies should be corrected for z (see e.g. Fenimore & 
Ramirez-Ruiz 2000). 
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correlation (within its scatter) for any z- As known GRB 980425 and GRB 031203 are also clear 
outliers of the Amati's relation (whereas GRB 030329 is consistent with it), i.e. of the strong 
associations 2/3 do not follow it (with respect to about 40 GRBs with known redshift). 

To conclude, within the limitations given by the (small) fraction of the GRB/SN sample 
we interestingly find that these events do not follow the Lag-Luminosity and the Peak Energy- 
Isotropic Energy correlations, in some cases for any assumed z. This - supported by the similarity 
with the behavior of GRB 980425 and GRB 031203 - appears to indicate that at least a few 
GRB/SN events form a separate class with respect to the majority of long GRBs. 

7. Summary and conclusions 

We presented a detailed analysis of the temporal and spectral properties of a sample of BATSE 
GRBs with an indication to be associated with a SN event, with the aim of determining whether 
these reveal any peculiar property with respect to the bulk of the detected GRBs. 

The "GRB/SN sample" was obtained combining GRB associations based on the presence 
of re-brightening in the afterglow light curves which might be interpreted as SN signatures and 
(more tentative) on positional/temporal coincidence from the cross correlation of GRB and SN 
catalogs. The latter ensemble was expanded with respect to previous works, by including the 
most recent and complete GRB and SN catalogs. The whole set comprises 36 GRBs (3 of which 
lack a corresponding SN). For GRBs with public data we performed a systematic analysis and 
compared their temporal and spectral properties with those of larger, "control" samples of BATSE 
bursts. 

The analysis of the temporal profiles revealed that in most cases the light curve presents a 
single peak. This property has been already considered by difl'erent authors (also on theoretical 
grounds) as a possible distinctive feature of GRBs associated with SN events (Bloom et al. 1998, 
Norris et al. 1999). The chance probability of finding the same fraction of single peaked bursts 
from a random extraction from the BATSE archive turned out to be ~ 0.4 %. 

We further analyzed the time profiles fitting a lognormal model and correlated the pulse 
parameters. Dim bursts turned out to have time lags larger than bright ones, as possibly explained 
by a larger viewing angle (loka & Nakamura 2001). This interpretation has also been invoked 
in the case of GRB 980425/SN 1998bw and GRB 031203/SN 20031w to account for their low 
luminosity compared to the bulk of GRBs (Waxman 2004). 

In terms of the spectral properties the GRB/SN sample shows softer spectra (both in terms 
of low energy spectral index and peak energy) with respect to typical bursts. This might be a 
distinctive feature of GRBs associated with SN as 2 out of 3 spectroscopic associations also 
share this atypical low peak energy. 

Finally, in the attempt to estimate their redshift, we considered two of the intrinsic correla- 
tions proposed for the population of long bursts, namely the Lag-Luminosity (Norris et al. 2000) 
and the £peak-£rad (Amati et al. 2002) relations. While some GRB/SN might be consistent with 
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them if not associated to the putative SN, some events do not obey such correlations for any 
redshift. This might indicate that they are effectively different or that these correlations only hold 
for particular sub-samples of GRBs (e.g. bright long bursts - see also Ghirlanda et al. 2004 who 
find that short bursts mostly populate a region above the iSpeak-^^rad correlation). In particular, the 
properties of this handful of GRB/SN again closely resemble those of 2 out of 3 spectroscopic 
associations which are also outliers of those relations. 

As discussed in the Introduction, previous results as well as the large positional uncertainties 
and wide time window for the selection of GRB/SN via catalog cross-correlation in fact point 
toward a low probability of finding true associations via this method. However, the above re- 
sults are suggestive of the opposite. Statistically, the GRB/SN sample contains (with respect to 
the whole of the BASTE events) an unusually high fraction of single peaked bursts, with softer 
spectra and some of those do not follow two of the proposed correlations for GRBs. Most of 
these properties are mimicking those of the spectroscopic GRB/SN associations. Although it is 
not obvious how to classify individual GRB/SN, in order to see any intersection among GRB/SN 
having 'unusual' characteristics, in Table 1 we schematically listed the three 'unusual' proper- 
ties for each GRB/SN (single-peakness, spectrum and consistency with the fipeak-^^rad relation). 
In 19 GRBs (other than GRB 980425) there are at least 2 properties which set the event as 
GRB 980425-like: 10 GRB/SN have 2 or 3 (3 of those are short and only 1 not selected by the 
catalog cross-correlation); 6 GRBs have 2 but not the third one (5 are short and 2 not selected by 
the catalog cross-correlation). 

A final indication that the catalog cross-correlations are not randomly sampling GRB and 
SN events comes from the properties of the selected SNe. We did not restrict our search for 
correlations on the basis of the SN Type, yielding a significant number of SNe la in the sample 
(see Table 1). While indeed it has been also proposed on theoretical grounds that also Type la SN 
might be associated to bursts (e.g. Berezinsky et al. 1996; Dar & De Rujula 2004), the current 
view (theoretically and observationally based) is that (some long) GRBs are likely associated 
with Type Ibc SNe. The intriguing finding is that among the GRB/SN events found by the catalog 
cross-correlations Type Ib,c SNe comprise about 28% of the events (i.e. 8/29), which is about 4 
times larger than expected from SN catalog statistic (see Valenti et al. 2005 who examined the 
Asiago SN catalog). This result also confirms the analysis independently carried out by Valenti 
et al. 2005 on BATSE catalog; indeed, these authors found an excess of Ib,c Types among SNe 
possibly associated with BATSE GRBs (~17%, with a significance level of ~ 93%). 

While the uncertainty on the reality of a single association remains high, these results indicate 
that: 1) there is a population significantly larger than that detected so far of GRB 980425-like 
events associated with SN; 2) the corresponding GRB have "atypical" temporal and spectral 
properties (and thus may suggest that not all GRBs are associated with at least typical 'SN'); 
3) they appear under-luminous in their observed y-ray emission (at the level of 0.1-1 %) with 
respect to "standard" GRBs. 
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All three aspects give statistical and physical information on the nature of GRB 980425-like 
events. The finding of a significant number of these under-luminous nearby bursts (e.g. viewed off 
axis) might in fact indicate from statistical arguments that these comprise the majority of GRBs 
(e.g. Sazonov et al. 2004). From the physical point of view, one promising scenario accounting 
for a qualitatively softer spectrum, a smooth, single-peak light curve and an apparently under- 
energetic event could be envisaged if the line of sight is at large angles with respect to the jet 
axis and the observed GRB emission is in fact prompt emission reprocessed (reflected) at such 
angles by circumburst material. A detailed study of the spectral and temporal properties expected 
in such case is ongoing. 

Estimates of the bolometric energetics of GRBs (e.g. Soderberg 2004) and the improved 
statistics which will be allowed in the HETE, INTEGRAL and starting Swift era offer the most 
obvious possibilities to test these results. 
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Table 1. The GRB/SN sample. The line separates GRBs with BATSE error boxes smaller and 
larger than 8 deg. (*) refers to a single-peaked GRB with substructure. There is no SN reported 
when there was only a SN bump present in the afterglow or in case of multiple (#) associations. 



® lAUC SN catalog. + Already in the Amati et al. (2002) sample; ' Wang & Wheeler (1998); ^ 
Hudec et al. (1999); ^ Piro et al. (1999); Germany et al. (2000); ^ Bloom et al. (1999); ^ Galama 
et al. (1998); ^ Holland et al. (2001); « Rigon et al. (2003); Terlevich et al. (1999). 



